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Abstract 

Crum's theorem and its modification a la Krein-Adler are formulated for the dis- 
crete quantum mechanics with real shifts, whose eigenfunctions consist of orthogonal 
polynomials of a discrete variable. The modification produces the associated polyno- 
mials with a finite number of degrees deleted. This in turn provides the well known 
Christoffel transformation for the dual orthogonal polynomials with the corresponding 
positions deleted. 

1 Introduction 

In a previous paper [1] we have developed a new paradigm for the orthogonal polynomials of 
a discrete variable [21 [Sj IH [5] , which is called the discrete quantum mechanics with real shifts. 
Various well-known orthogonal polynomials, for example, the (g)-Racah, (g)-(dual)-Hahn, 
(g)-Meixner, etc. [5] are derived in p] as the main parts of the eigenvectors of specific tri- 
diagonal real symmetric (Jacobi) matrices of finite or infinite dimensions. The orthogonality 
of the eigenpolynomials is evident by construction and the orthogonality weight function is 
obtained as the square of the groundstate eigenvector. The dual pair of orthogonal poly- 
nomials {the Leonard pair) [6l [7] is naturally defined and the three term recurrence relation 
for the dual polynomials is built in automatically. It should be emphasised that the present 
paradigm offers a unified understanding [HI |9] of the (g)-Askey scheme of hypergeometric 
orthogonal polynomials [21 [5] . Various properties of the discrete quantum mechanical sys- 
tems with real shifts, for example, the duality, the orthogonality, the exact solvability in the 



Schrodinger as well as the Heisenberg picture [10], the dynamical symmetry algebras (the 
Askey- Wilson algebra |Tl] and its degenerations [I2]), the quasi-exact solvability, etc. can 
be understood in a unified fashion, not one by one specifically. 

In the present paper we start with the formulation of the discrete quantum mechanics 
version of the well-known theorem by Crum [13] for one-dimensional quantum mechanics or 
Sturm-Liouville theory. Crum's seminal paper asserts, in the language of quantum mechan- 
ics, the existence of an infinite family of associated Hamiltonian systems {i.e. the Hamilto- 
nian together with its eigenvalues and eigenfunctions) which are essentially iso-spectral with 
each other. We will show that the situation is the same in the discrete quantum mechan- 
ics with real shifts. The associated Hamiltonians are tri-diagonal real symmetric (Jacobi) 
matrices, too. See Fig. 1 for the general structure of the associated Hamiltonian systems. 
Crum's theorem for the discrete quantum mechanics with pure imaginary shifts [Ul [151 E] 
has been formulated in [16]. This covers the orthogonal polynomials with absolutely contin- 
uous weight functions, for example, the Wilson and the Askey- Wilson polynomials and their 
various degenerate forms. See also [17] in this context. 

Next we will present the discrete quantum mechanics (with real shifts) version of Krein- 
Adler [18] modification of Crum's theorem [19]. The Krein-Adler modification goes as fol- 
lows: starting from a given Hamiltonian with eigenvalues {S{n)} and the corresponding 
eigenfunctions n = 0, 1, . . ., one can again construct an associated Hamiltonian system 

by deleting a finite number of energy levels {S{dj)}, dj G Z>o, j = !,...,£, which satisfy 
certain conditions (14. ip . We will demonstrate that the situation is the same in the discrete 
quantum mechanics. The associated Hamiltonian is again a tri-diagonal real symmetric (Ja- 
cobi) matrix and the eigenfunctions form a complete set of orthogonal functions. See Fig. 2a 
for its schematic structure to be contrasted with the original Crum's case in Fig. 2b, which 
corresponds to the very special choice of deletion {di, d2, . . . , di} = {0,1,2, i — 1}. If the 
starting system is exactly solvable, as those examples given in [T], the modified system is 
also exactly solvable. 

Crum's theorem and its modification will be presented generically, at the beginning. 
That is, no extra condition is imposed on the functions B{x) and D{x) in the Hamiltonian, 
other than the positivity, the boundary (asymptotic) conditions (12.31) . Many formulas will 
be drastically simplified when the eigenfunctions consist of polynomials and/or the system 
is shape invariant [2U], which is the case for most practical applications. We will present 
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these simplified formulas, too. Due to the lack of the generic 'oscillation theorem' in discrete 
quantum mechanics, we do not have categorical proof of the hermiticity of the Hamiltonian 
% (14.391) obtained by the application of the modified version of Crum's theorem. For each 
particular case of specific polynomials the hermiticity is verified. 

It is important to stress that the modified Crum's theorem provides a unified theory 
of dual Christoffel transformations [21j for orthogonal polynomials of a discrete variable. 
The dual Christoffel transformation has the merits that the formulas are universal, concise 
and algorithmic compared with the original Christoffel transformation, which should be 
performed specifically for each case, the polynomials and the set of deletion D. This situation 
is explained in some detail in section [61 

As an illustration of the modified Crum's theorem, we will derive in Appendix uni- 
fied expressions of the eigenpolynomials for the special case of deletion {di, (i2, . . . , rf^} = 
{1, 2, . . . , £}. The same formulation and examples for the discrete quantum mechanics with 
pure imaginary shifts are reported in [19] . As is well known these orthogonal polynomials of 
a discrete variable have many important applications in many arenas of physics and mathe- 
matics OEIH]. To name a few recent applications, the linear quantum registers [22j, quan- 
tum communications [23| and the birth and death processes. As shown in [23], the explicit 
examples of 18 orthogonal polynomials in [1], the (g)-Racah, (g)-(dual)-Hahn etc, provide 
exactly solvable birth and death processes [H [25]. That is, for the given birth and death 
rates {B{x),D{x)} which define the Hamiltonian (12. ip . the corresponding transition proba- 
bilities are given explicitly, not in a general spectral representation form of Karlin-McGregor 
[2S] • By applying the present modification of Crum's theorem to these polynomials, one can 
generate an (in)finite variety of exactly solvable birth and death processes. 

This article is organised as follows. In section [2] we recapitulate the essence of the discrete 
quantum mechanics with real shifts in order to introduce necessary notions and notation. 
In section [31 the discrete quantum mechanics version of Crum's theorem is formulated in 
its full generality. The modification of Crum's theorem a la Krein-Adler is developed in 
section [H Section [S] provides various simplifications of the formulas in the presence of shape 
invariance and/or polynomial eigenfunctions. In section [SI the Christoffel transformations 
for orthogonal polynomials of a discrete variable in discrete quantum mechanics are shown 
to be dual to the modification of Crum's theorem developed in the preceding sections. The 
final section is for a summary and comments. Appendix gives the simplest examples of 
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the modified Hamiltonian systems obtained by deleting the lowest lying ^ excited states for 
various exactly solvable Hamiltonians discussed in [1]. 

2 Discrete Quantum Mechanics with Real Shifts 

Let us review the discrete quantum mechanics with real shifts formulated in [1] . The Hamil- 
tonian % = {'Hx,y) is a very special type of tri-diagonal real symmetric (Jacobi) matrices and 
its rows and columns are indexed by non-negative integers x and y {x,y = 0, 1, . . . ,Xmax), 
in which Xmax is either finite (xmax = N) or infinite (xmax = oo). The Hamiltonian "H has a 
form 

n = -^/B{^e^^/D{^- ^/D{^e-^^/B{^ + B^^^ (2.1) 
n.,y = -^/B{x)D{x+l)S,+,,y - ^B{x-l)D{x)5,^i^y + {B{x) + D{x))5^^y, (2.2) 

in which 9 = ^ is the differentiation operator (df){x) = {e^^f){x) = f{x ± 1). The 

two functions B{x) and D{x) are real and positive but vanish at the boundary: 

B{x) > 0, D{x) > 0, D{0) = ; B{x^^^) = for the finite case. (2.3) 

The problem at hand is to find the complete set of eigenvalues and the corresponding eigen- 
vectors of the hermitian matrix H (n^ax = N or oo): 

n(j)n{x) = £{n)(j)n{x) (n = 0, 1, . . . ,nmax), (2.4) 

which is the Schrodinger equation of the discrete quantum mechanics with real shifts. It has 
non-degenerate spectrum thanks to the generic property of Jacobi matrices, 

£{0) <£{1) <S{2) < ■■■ . (2.5) 

The Hamiltonian (12. ip can be expressed as a product of a bi-diagonal lower triangular matrix 
and a bi-diagonal upper triangular matrix A: 

T~L A i^Ax^y)^ Al^ ((^^)a;,?/) i^Ay^x) (-^j V 0, 1, ... , Xjnax); (2'6) 

A= ^/B(^-e^^/D{x), A^ = ^B{x) - ^D{x) , (2.7) 

Ax,y = y^B{x) 6x,y - y^D{x+l)6x+l,y, {A^)x,y = \/B{x) 6x,y - ^/D{X) 6x^l,y. (2.8) 
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The zero mode Acpoi^x) = is easily obtained: 

00 (a;) 



with the normahsation 0o(O) = 1 (convention: nfc=n * = !)■ Note that the groundstate wave- 
function (po{x) of (12. ip is positive {(Pq{x) > 0) throughout the range of x (x = 0, 1, . . . , Xmax)- 
Then the Hamihonian is positive semi- definite: 

^(0) = 0. (2.10) 

The finite £^-norm condition of (poi^) 

Y,M^Y<oo, (2.11) 

imposes conditions on the asymptotic forms of B{x) and D{x) for the infinite case. 
The orthogonahty relation is 

2^ max 

(0n,0m) = ^0n(x)0m(x) = (u, m = 0, 1, . . . , Wmax) • (2.12) 

Here l/ci^ is the normalisation constant. For the infinite ), the eigenfunctions 

should satisfy the asymptotic condition 0„(xmax) = and the finite £^-norm condition. 
If the eigenfunction has the factorised form, 

M^)=M^)Pnivi^)), (2.13) 

where Pniv) is a polynomial of degree n in the sinusoidal coordinate rj = rj{x) with the 
normalisation, 

P„(0) = 1 (n = 0,l,...,n„,ax), P--i(r/)=0, 7^(0) = 0. (2.14) 

The sinusoidal coordinate r]{x) is the central notion of exactly solvable quantum mechanics 
[HEmn]. it undergoes sinusoidal motion through the time evolution governed by the Hamil- 
tonian "H, see (2.28) of [TU]. It is also established that the Heisenberg operator equation for 
the sinusoidal coordinate rj{x) is exactly solvable PQ[Hl[in]- In the discrete quantum mechan- 
ics with real shifts, there are five sinusoidal coordinates, see (4.72)-(4.76) of [Tj. Then the 
above orthogonality relation becomes that of the orthogonal polynomials 

^max 

^ 0o(a;)^P„(?7(a;))Pm(?7(x)) = — (n, m = 0, 1, . . . , n^ax)- (2.15) 
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Here the discrete weightfunction (poix)"^ is given explicitly (12. 9p . The orthogonal polynomials 
Pn{ri{x)) are the 'eigenfunctions' of the similarity transformed Hamiltonian 7/ = {'Hx,y) 

(x, y 0,1,..., Xjnax) 

n oHo 0o(x) = B{x){l - e^) + D{x){l - e"^), (2.16) 

T-L^^y = B{x){6cc,y - Sx+i,y) + D{x){6r,^y - (2.17) 
nPM^)) = Sin)PM^)) (n = 0,l,...,n^ax) (2.18) 

= B{X) {Pnivix)) - Pnivix + 1))) + D{X) (Pnivix)) ' Pnivi^ " 1))) • 

The dual polynomial Qx{£), which is a degree x polynomial in S with the normalisation 
condition 

Q,(0) = 1 (x = 0,l,...,w), Q-ii£)=0, £iO) = 0, (2.19) 
is defined by the three term recurrence relation 

B{x){Q,{£) - Q.+ii£)) + D{x) {Q,{£) - Q.-i{£)) = £Q,{£) (x = 0, 1, . . . , w). (2.20) 

For £ = £{n), the above difference equation for Pn{ri{x)) (I2.18P is identical with the three 
term recurrence relation under the identification 

Pnivix)) = Qxi£in)) (X = 0, 1, . . . ,Xmax), (^ = 0, 1, . . . , n^ax) • (2.21) 

This establishes the duality [6| [?] . The orthogonality relation for the dual orthogonal poly- 
nomials Qx{£{n)) takes a dual form to (12.151) 

J2 dlQxi£in))Qyi£in)) = —— b^y {x,y = 0, 1, . . . ,x^ax). (2.22) 

For more details of the duality in the context of discrete quantum mechanics with real shifts, 
see § 3 of [1] . 

For all the examples presented in [JLJ, the eigenfunctions have the form (I2.13P and the 
square of the groundstate wavefunction (pQ^x)"^ can be analytically continued to the whole 
complex X plane and it has (at least) a simple zero at integral points outside [0,Xmax], 

0o(x)' = O forxGZ\{0,l,...,w}, (2.23) 

due to the factors of (poi^xY such as (g ; q)~^ = {q^~^^; q)oo{<l ', q)^, <i)x, etc. Thanks to 

this situation, various expressions in the subsequent sections whose arguments appear to go 
beyond the defined range of [0, Xmax] cause no harm. 
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3 Crum's Theorem 



Crum's theorem [13] describes the relationship between the original and the associated Hamil- 
tonian systems, which are iso-spectral except for the lowest energy state. The relationship 
among various associated Hamiltonian systems is quite general, as depicted in Fig. 1, and it 
is shared by the most general ordinary quantum mechanical systems [13] as well as discrete 
quantum mechanical systems with real and pure imaginary shifts [T6l [T7]. As shown below, 
the factorised form of the Hamiltonian fl2.6p is essential. 

3.1 Deletion of the ground state 

For later convenience, let us attach the superscript ^ to all the quantities of the original 
Hamiltonian system, = = = A, 4^f ^t^ b'^'^^x) = B{x), 

D^^^{x) = D(x). Let us define an associated Hamiltonian T/'^^ by simply changing the order 
of A and A^: 

^[i]drf^[o]^[o]t_ (3_1) 
By construction A and A^ intertwine Ji^^^ and H'^^: 

^[o]^[o] ^ ^[i]^[o]^ ^[o]t^[i] ^ ^[o]^[o]t^ (3^2) 

The matrix elements of the associated Hamiltonian "Ht^l are 

= -^/B{x + l)Dix + 1) - ^/Bix)Dix) + D(x + l)(l -(5,,,_J)4,y 

(3.3) 

For the finite case (xmLx = N, B{N) = 0), the matrix elements in the (A^ + l)-st entry 
vanish, namely 'H'^l has the following form: 



(3.4) 



where H^^'' is an x matrix and is an A^- dimensional zero column vector. Therefore the 
eigenvalue problem of T-L^^^ reduces to that of T-L^^^' and the trivial one dimensional part with 
zero eigenvalue corresponding to the deleted groundstate of the original Hamiltonian. To get 
rid of this trivial part we define Xmlx as Xmlx = A^ — 1. For the infinite case (xmax = oo), Xmlx 
is defined as Xmax = oo. To treat the finite and the infinite cases in parallel and to avoid 
notational cumbersomeness, we write T-L^^^' as "H'^^ by adopting the following convention. 
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Namely H'^^ represents an (xmax+ 1) x (a;max+ 1) matrix H^^^ = (hL y) (a;, y = 0, 1, . . . , XmL) 
with XmLx = — 1 or oo. 

We will show that the associated Hamiltonian system H^^^ is iso-spectral to the original 
Hamiltonian system 'H^'^^ and the eigenfunctions are in one to one correspondence, except for 
the groundstate. Thanks to the first of the above relation fl3.2p . it is trivial to verify that 
the eigenfunctions of the associated Hamiltonian system "H'^' are generated algebraically by 
multiplying to the eigenfunction of the original system: 

= (n = 1,2, . . . ,n^,x), (3.5) 

n^'^(P\l\x) = S{n)<f)^^\x) {n = l,2,...,n^^). (3.6) 

The orthogonality relation is 

■''max 

(0™, 01^1) = E 0L''(^)0S(^) (n, m = 1, 2, . . . , n_) 

x=0 

= (^M0L°U[°i0[;ii) = (0[°ut°iUM0!°^) = (0L°U[°i0l^i) = ^(n) (3.7) 

For the finite case this gives all the eigenvectors of H^^^. For the infinite case, suppose 
the associated Hamiltonian H^^^ has an eigenfunction (p'i^x) with an eigenvalue £' other than 
those listed above: 

?^W0'(x) = S'(j)'{x). (3.8) 
Again, thanks to the second of the relation (13. 2p . it is trivial to verify 

?{[°U[°]V'(a;) = ^U[°]V'(a;). (3.9) 

Due to the completeness of the spectrum of the original Hamiltonian 7/'°], the provisional 
eigenvalue S' must belong to the spectrum (12. 5p S{n) for n = 1,2. . . . , nmax- In other words, 
S' cannot be vanishing, S' ^ 0. Suppose that is the case {S' = 0), then 0' is annihilated 
by A'^^^^. But it is easy to see that there exists no non-zero solution (the finite case) or 
no finite £^-norm solution (the infinite case) of the equation A^'^^^ (f)' (x) = 0. Thus we have 
established that the associated Hamiltonian system "H^^' is essentially iso-spectral to the 
original Hamiltonian system 'H^'^^ and the eigenfunctions are in one to one correspondence, 
except for the groundstate with the wavefunction 0q°'(x). 
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If the groundstate energy £^(1) is subtracted from the associated Hamiltonian "Hl^l, it is 
again positive semi-definite and can be factorised as above: 

= ^WUW+^(1), (3.10) 
A^^ = ^JB\%x)-e^^JD[%x), A^^^ = ^J B^^x) - sj D[%x) e'^ , (3.11) 



B^\x) \Ib\%x + l)D[o](x + 1) ^lljj±l) , (3.12) 

01 (^) 



D^'\x) = Ji?M(x)DM(x) '^^ (3.13) 



def /„rnl/ X r^rnw n 4>l\x - 1) 

0?'(^) 

Since 0^' is the groundstate of 'H^^\ it does not vanish inside the interval [0,Xmax]. Thus 
B^^^(x) and D^^^(x) are non-singular and positive. Note that D^^^(0) is not yet defined because 
(f)f\—l) is not defined. Due to the factor D^^^{x) in (13.131) and property D^^^{0) = 0, we define 
/}[i](0) =^ 0. In concrete examples the expression of (/)f'(x)^ can be analytically continued and 
it vanishes at x = —1. This fact also supports the definition D^^^(0) = 0. The groundstate 
wavefunction (pf^ (x) has the following form 



0W(a;) = ^M0fl(x) = ^5M(x)0fl(x) - y/DM(x + 1) 0fl(x + 1)(1 - 5^^jo,J. (3.14) 
Due to the vanishing factor B^^^{xma.x) = for the finite case and (f)n\xma.x) = for the 



infinite case, this expression of the groundstate wavefunction 0f'(a;) vanishes at a; = 



max. 



0f'(i) = O. (3.15) 

Thus we have B^^\xma.x) = for the finite case. 

It is easy to verify that 0^^' (x) is annihilated by A^^^ , 

^W0W(a;) = O. (3.16) 

The equation (13.101) is shown by elementary calculation. 

3.2 Repetition 

Starting from (13.101) . the second associated Hamiltonian system 'H^^l can be defined by re- 
versing the order of ^'^^^ and A^^^. This process can go on repeatedly. 
Here we list the definition of the s-th quantities step by step for s > 1, 

^Wdcf^[.-i]^[.-i]t + ^(,_l)^ xl^L = iV-sor oo, (3.17) 
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"=^' {n = s,s + l,..., n^ax), (3.18) 

= ^5W(a;) - e^^Jol^), = ^5^(3.) _ ^D[-\x)e~^, (3.19) 

5H(x) JE[--il(a: + l)Di^-\x + 1) '^''^['^ ^ , (3.20) 

0r (x) 

dW(x) ^=yi?[-^(x)D[-il(x) '^^^'il'^^ (3.21) 

0s (x) 

Like in the s = 1 case we set DW(0) = 0. Recall that is an (Xmax 

+ 1) X (xl^L + 1) 

matrix "Ht*' = ('HL'^,'j/) (x, ?/ = 0, 1, . . . , Xmax). Then we can show the following for s > 0, 

?^W</.W(x)=£:(n)0H(x) (n = s,s + l,...,n^ax), (3.22) 

^W0H(x) = O, (3.23) 

-HH =^WUH+£(s), (3.24) 



S — 1 



E '^n^(^)'^m (^) = n^^H - ^(^■)) ■ (n, m = s + 1, . . . , n„ 



x=o i=o " 

(3.25) 

and -B'*'(a;max) = for the finite case. We have also for s >1 

<Pt\^) = \ ^( TT'^«(^) (« = s,s + l,---,^max). (3.26) 

t[n) — t[s — 1) 

The situation of Crum's theorem is illustrated in Fig. 1. 

As in the original Crum's case [13], the eigenfunction (j)n\x) can be expressed in terms 

of determinants. Let us define the Casorati determinant for n functions fj{x) as 

W[/i, . . . , U]{x) = det (Mx + J - 1)) , (3.27) 

(for n = 0, we set W[-](x) = 1), which satisfies 

n-l 

W[gh, gf2, gfn]{x) = l[g{x + k)- W[/i, /a, . . . , /„](x), (3.28) 

fc=0 

W[W[/i,/2,...,/n,^7],W[/i,/2,...,/„,/l]](x) 

= W[/i, /2, . . . , fn]ix + 1) W[/i, /2, . . . , fn, g, h]{x) (u > 0). (3.29) 
By using the Casorati determinant, we obtain 

1=0 W[0o,01,...,0s-l](x + l) ^ ' 
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n^i] n^2] ^[3] 

Figure 1: Schematic picture of Crum's theorem 

4 Adler's Modification of Crum's Theorem 

Crum's theorem describes the construction of an associated Hamiltonian system which is 
iso-spectral to the original one with the lowest energy state deleted. Adler's modification 
[18] of Crum's theorem is the construction of an associated Hamiltonian system which is 
iso-spectral to the original one with finitely many states deleted. 

Let us choose a set of i distinct non-negative integers l] V =^ {di,d2, ■ ■ ■ ,di} C Z>o, 
satisfying the condition [18] 



JJ(m - dj) > 0, Vm G Z>o. 



(4.1) 



^ Although this notation dj conflicts with the notation of the normalisation constant dn in p.l2p . we 
think this does not cause any confusion because the latter appears as ^ Smn- 
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This condition means that the set T> consists of several clusters, each containing an even 
number of contiguous integers 

dki-, c^fci + 1, ■ ■ ■ ) ; c^fes, c^fca + 1, ■ ■ ■ ) j c^fcs; dk^ + 1^ - ■ ■ ^dk^ ] ■ ■ ■ , (4.2) 

where ci^j + 1 < rf^g, dk^ + I < dk^, ■ ■ ■ . If dk^ = for the lowest lying cluster, it could 
contain an even or odd number of contiguous integers. The set V specifies the energy levels 
to be deleted. Deleting an arbitrary number of contiguous energy levels starting from the 
groundstate {V = {0,1,2, i — 1}) is achieved by Crum's theorem discussed in §[3J 



da 
d-> 



d^ 
ds 

di 



it 



A 





original 

system 



* 

modified 
system 



A 



2 
1 


original 

system 





after i steps 
of Crum 



^di...di ' ' ' A(iid2Adi 

Figure 2a: Generic case 



Figure 2b: Crum's case 



We will construct associated Hamiltonian systems corresponding to the successive dele- 
tions Tidi... (and Ad^..., A]i_^ , etc.) step by step, algebraically. It should be noted that some 
Hamiltonians in the intermediate steps could be non-hermitian. 

4.1 First step 

For given di, the original Hamiltonian 7i can be expressed in two different ways: 



n = AU = AlAd,+S{di), 

Ad/= y^BZi^-e'y^DZi^, AI/= ^/B^i^) - y^D^i^ e'^ 



(4.3) 
(4.4) 
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BM'=' ./mWTT)^^^^^^, DM'^^B{x-l)D{x) '^'f. (4.5) 
and we have 

AdMx) = 0. (4.6) 

As in §[3]we set Dd^ (0) = 0. Note that Bd^ (a;max) = for the finite case. Unless di = 0, Bd^ (x) 
and Dd^ (x) are not always positive due to the zeros of (pdi (x) ■ It is important to note that 
•^d-^ in (14.41) is a 'formal adjoint' of Ad^ due to the above mentioned sign changes of Bd^x) 
and Dd^x). We stick to this notation, since the algebraic structure of various expressions 
appearing in the deletion processes are best described by using the 'formal adjoint'. By 
changing the order of Ad^ and , let us define a new Hamiltonian system 

Hd.'^AdAi+^id,), xtJ=N-lov^. (4.7) 

As in §131 the Hamiltonian "H^^ represents an ix'^^-^ + 1) x (x^^^^ + 1) matrix "H^^ = {'Hdi;x,y) 
{x,y = 0,1, . . . , x^a^-^). It is easy to show that the 'eigenfunctions' of this Hamiltonian are 
given by 

4>dinix)= Ad^4>n{x) (nG {0, l,...,nmax}\{'^l}), (4.8) 

T-id^4>d^n{x) = £{n)(j)din{x) (n G {0, 1, . . . ,nmax}\{c?i})- (4.9) 

Thus the energy level di is now deleted, (t)d]_dAx) = 0, from the set of 'eigenfunctions' 
{0(iin(x)} of the new Hamiltonian "H^i- 

4.2 Repetition 

Suppose we have determined the Hamiltonian 'Hdi...ds together with the eigenfunctions 
(pdi ...dsn{x) with s deletions. They have the following properties: 

nd,...d. =Ad,...dAi,„d^+£ids), xt-'^ ""^'N-s or oo, (4.10) 
Ad^ ...d, = ^Bd,...ds{x) - e^^Dd,...ds{x), 

= VBd....dAx)- VDd,...dMe-', (4.11) 

jBd,...d..Ax+l)Dd,...d.^Ax + l) '^l'-'^^^!/^ (^>2) 
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0di...ci.n(2;) = ^di...(i>di...<i._in(a;), (4-14) 

'Hdi...ds<l>di...d,n{x) ^ £{n)(f)di,„d,n{x), (4.15) 

where n G {0, 1, . . . , ninax}\{o?i, ■ ■ ■ , o^s}- As before we set -Ddi ...(i^(0) — 0. We also have 

Bdi ...(i„(a^max) = for the finite case. We note that the following relations hold: 

Bd, ...dXx)^Dd,... d. (x + 1) ( '^'^.-'^ + ) ' (.>1), (4.16) 

R . ^n ( [ Bd,...d._,{x^\)Dd,...d,_,{x^\) {s>2) 

5.......(x)i?.......(x + l) = |^^^^^^^^^^ (. = 1) ' ^^-1^) 



^di ... (i« (2^) + ^di ... (i« (a;) + ^ (4) 



We have also 



...d._,{x) + L><i, ...d,_,(x + 1) + £:(d,_i) (s > 2) 

S(x) + D{x) (s = 1) ' 

(4.18) 



^d^...d,-^n{x) = Z'"'^^; , ^ 0di...d,n(a;) (n e {0, l,...,nmc«}\{c?i,...,(ij). (4.19) 

Next we will define a new Hamiltonian system with one more deletion of the level o?s+i. 
We can show the following: 

y-di-ds = '^di...dsds+i'^di...d,ds+i +£{ds+i), Adi...dsds+Adi...d,d,+i{x) = 0, (4.20) 

Adi...dsds+i = \/ Bdi...dsds+i{x) - ^^\J Ddi...dsd,+iix), 

= ^/Bd,...d. d.^, (x) -^Dd,...d. d.^, (x) e-^ (4.21) 

Bd,M,A^) - ^jBd,...dA^+l)Dd,...dA^ + ^) '^Y'''"'^'^^^^ ^ (4-22) 

(Pd^...dsds+AX) 



<Pdi...dsds+i{x - 1) 



Dd....d.d.,.{x) VBd....dA^)Dd,...dM '^"'Y'"''''^ (4.23) 

(Pdi...ds d,+Ax> 

These determine a new Hamiltonian system with s + 1 deletions: 

= ^d,...d,+,4,...,,^, +^(4+i), xt^'^+' 1 or 00, (4.24) 

(Pdi...d,+in{x) '^^^ Ad^...d,+Adi...dsnix), (4.25) 

Ud^...ds+i(t>di...ds+in{x) ^ £{n)(f)d,...d,+inix), (4.26) 
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where n e {0, 1, . . . , n^s.^}\{di, 

By induction we can show that the eigenfunction is expressed in terms of the Casorati 
determinant: 



.,,„(a;) = (-ir 



^di ...d, X + s + 1 - A;) • — — ^- , u ^ ■ 



(4.27) 
(4.28) 



We can also show the following: 

k=l 

s 

YlDa,...ci,{x + s + l-k) 



flB{x + k- l)D{x + k) . (4.29) 

\ k=l 



k=l 



W[0,„...,0,J(x) 
W[0di,...,0dj(a;) 



, Y\B(x + k-l)D{x + k) ^^^.^ ^ 
\t}i W[(/.,„...>dJ(^ + l) 



. (4.30) 



4.3 Last step 



After deleting all the T> — {di, • • • , dg] energy levels, the resulting Hamiltonian system 
?^x) = y-di...de, -A-v = Adi...de^ etc has the following form: 



Hv = AvAl, + S{de) , x^^ = TV - ^ or oo, 

f dcf 



X> def 



Bv(x) - J DM e-^ 



Av = Vb^ - e^/D^, A 
Bv{x) = ySdi...d,_i(2; + l)£>di...d,_i(2; + 1) — ^^^^ 



def 



- 1) 



(pVnix) = ^X'0di...d,_in(a^) (JT- £ {0, 1, . . . , n^axlX^) , 
'Hv<i>vnix) = £{n)(t)vn{x) (ne {0, 1, . . . , nmax}\^^)- 



Now 1-Lx) has the lowest energy level ji: 



ji =^ min{n | n G {0, 1, . . . , n^s^}\D}, 



(4.31) 
(4.32) 

(4.33) 

(4.34) 

(4.35) 
(4.36) 

(4.37) 



with the groundstate wavefunction (p^^x) 



4>i^{x) 4>Vti{x) = 4>di-deti{x). 



(4.38) 
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Then the Hamiltonian system can be expressed simply in terms of the groundstate wave- 
function which we will denote by new symbols 7{, A, etc: 



■H = 'Hv = A^A + S{n), 

def 



- £ or cx), = 0, 



A ^ Av, = ^B{x) - e'^Dix), A^ ^ A^^ B{x) 
B{x) = Bv^{x) = y^Bj,{x + l)Dj,{x + 1) "^i"" ■ 



D{x) e 



D{x) = Dv^,{x) = ^Bv{x)Dv{x) 



1-L(t)n{,x) = 8{n)(j)n{x), (f)n{x) = (f)vn{x) {u E {0, 1 , . . . , n^i,^}\V) , 
B{x) + D{x) + Sifi) = Btj{x) + Dv{x + 1) + £{di). 



(4.39) 
(4.40) 

(4.41) 

(4.42) 

(4.43) 
(4.44) 



Again we set -D(O) = 0. We have -B(xmax) = for the finite case. Note that the final 
Hamiltonian 1-L is an (a^max + 1) x (xmax + 1) matrix l-i = {T-Lx,y) (a;, y = 0, 1, . . . , Xmax)- 

In terms of the Casorati determinants we obtain the expressions of the final eigenfunction 
(j)n{x) and the final functions 13{x) and D{x) {I > 0) 



= (-1/ 



p . N W[0rf,, . . . ,(j)d^,(j)n]{x) 

'^•••'^^'^^■W[0,„...,<^.j(x + 1) 



\ k=l 



di ...d^ 



{x + i + l-k)- 



B{x) 
D{x) 



^B{x + e)D{x + £ + l) 



W[0d,,...,0dj(x) W[0d,,...,0d,,0^.](a; + 1) 



'B{x-l)D{x) 



W[(f)d„ 0dj(x + 1) W[0d,, . . . , 0^](a;) 

W[0d,,...,0rfj(x+ 1) W[0rf,,...,0rf„0^](x- 1) 



W[0rf,,...,(/)rfj(x) W[0rf,,...,0rf,,0^](x) 

From f lCTjl -f iOOj) we also have {£ > 0) 



(4.45) 

(4.46) 

(4.47) 
(4.48) 



YlBd,...d,{x) 

k=l 

i 

l[Dd,...d,{x + e + i-k) 



. YlB{x + k-l)D{x + k) 

\ k=l 



k=l 



. YlB{x + k-l)D{x + k) 

\ k=l 



W[0rf,,...,0rfj(x + 1) 

W[0d,,...,0dj(x) 
W[0rf,,...,0dj(x) 

W[0rf,,...,0dj(x + 1)' 



(4.49) 



(4.50) 



Therefore the functions B{x), D{x) and 0„.(x) are symmetric with respect to c?i, . . . , di, and 
the final Hamiltonian "H is independent of the order of {dj}. 
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Let us state the discrete quantum mechanics version of Adler's theorem; // the set of 
deleted energy levels V = {di, . . . ,di} satisfy the condition (14. ip . the modified Hamiltonian 
is given by'H = Tidi ...di = A'^A+Sl^) with the potential functions given by (I4.47P - fl4.48l) and 
its eigenf unctions are given by fl4.45p - fl4.46p . The discrete QM version of Crum's theorem 
in §[3] corresponds to the choice {di, . . . , di} = {0,1, ...,£— 1} and the new groundstate is 
at the level fi = i and there is no vacant energy level above that. We have not yet proven 
the hermiticity of the resulting Hamiltonian "H and the reality of the eigenfunctions (f)n{x) 
categorically for the discrete quantum mechanics, even when the condition (14. ip is satisfied 
by the deleted levels. This is due to the lack of the difference equation analogue of the 
oscillation theorem. It should be stressed that in most practical cases, in particular, in the 
cases of polynomial eigenfunctions, the hermiticity of the Hamiltonian is satisfied. 

The orthogonality relation of the complete set of eigenfunctions is 

(0n,0m) = ^4>n{x)(j)m{x) = Y\_{^i^) - ^{dj)) ■ ^nm {n,m G {0, 1 , . . . , nmax}\^^) ■ 

(4.51) 

Note that the coefficient of 6nm is positive for V satisfying (14. ip . 

5 Simplifications 

In the preceding sections, Crum's theorem and its modification are presented in their most 
generic forms. These formulas are simplified substantially when the system is shape invariant 
[2U] and/or the eigenfunctions consist of polynomials. 

5.1 Formulas for Crum's theorem 
5.1.1 Shape invariance 

Shape invariance is a well known sufficient condition for exact solvability in ordinary quan- 
tum mechanics. The situation is exactly the same in discrete quantum mechanics. Shape 
invariance simply means that the s-th associated Hamiltonian system has the same form as 
the original one with certain shifts of the contained parameters. 

The Hamiltonian may contain several parameters A = (Ai, A2, . . .) and we indicate them 
symbolically as "H = 'H(A), A = ^(A), £{n) = £{n;X), (pn{x) = 0„(x;A), etc. Let us 
consider the case that the potential functions of the first associated Hamiltonian B^^^(x) = 
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B^^^{x ; A) and D^^^(x) = D^^^(x ; A) have the same forms as the original functions B and D 
with a different set of parameters and up to a muhiphcative positive constant k G ]R>o: 

A) = A'), D^^\x ; X) = kD{x ; X'). (5.1) 

Here the new set of parameters A' is uniquely determined by A. In concrete examples we 
can choose an appropriate set of parameters such that A' = A + 5, with a shift of parameters 
denoted by 6. In the following we assume this. In fact, (15. ip is equivalent to the definition 
of the shape invariance by Odake and Sasaki [HI [H [H] , 

^(A)^(A)^ = kA{X + 5)U(A + 6)+ £{!■ A). (5.2) 

It should be stressed that the above definition of shape invariance is more stringent and 
thus more constraining than the original definition of Gendenshtein [20]. This condition is 
rewritten as 

^/B{x + 1; X)D{x + 1; A) = B{x] A + 6)D{x + 1; A + 6), (5.3) 
B{x] A) + D{x + 1; A) = k{B{x] X + 6)+ D{x; X + d))+ £{1; A). (5.4) 

The shape invariance condition (15.21) combined with Crum's theorem implies 

A^'\X) = K^A{X + sS), ^W^(A) = fi:M(A + s(5)^ (5.5) 
n^'\X) = k"H{X + sS)+£{s]X), (5.6) 

and the entire energy spectrum and the excited states wavefunctions are expressed in terms 
of £^(1; A) and (t>o{x ; A) as follows: 

n-l 

£(n;A) = ^/€^£(l;A + s5), (5.7) 

s=0 

0„(x ; A) oc ^(A)U(A + 5)U(A + 26^ • ■ ■ ^(A + (n - 1)5) Vo(a; ;X + nS). (5.8) 

Therefore the shape invariance is a sufficient condition for exact solvability. We have also 

A{X)M^; A) = /„(A)0„„i {x; X + 6), (5.9) 

AiXy^n-i{x;X + d) = /B(O^6„_i(A)0„(x;A), (5.10) 

where /„(A) and 6„_i(A) are the factors of the energy eigenvalue, S{n; A) = /„(A)6„_i(A). 
It is interesting to note that the polynomial eigenfunctions are not the direct consequence 
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of the shape invariance. When the eigenf unctions consist of polynomials (I2.13p . the above 
formula (15.81) could be called the Rodrigues formula for the polynomials. Relations (15. 9p - 
(I5.10p would translate into the forward and backward shift relations [H |8] for the polynomial 
eigenfunctions. 

5.1.2 Polynomial eigenfunctions 

Here we consider a generic Hamiltonian (12. ip . That is the shape invariance is not assumed. 
Let us define a function ri{x) as a ratio of 0i(x) and (poi^), 

^ = a + br^ix), (5.11) 

(Pq{x) 

where a and b {b 0) are real constants. Although rj{x) is not well defined without specifying 
a and b, this ambiguity (affine transformation of ri{x)) does not affect the following discussion. 
Then fl312D and fl3J3D imply 

Bl'lM = B(. + l) '''" + ^7'" + '' . jt'lM^CM "'"-"-"'"; . (5.12) 

Let us assume further that the n-th eigenfunction (pn/'Po is a degree n polynomial in this 
ri{x) for all n > 2: 

n 

0„(a;) = 0o(a;)P„(?7(a;)), P„(?/) = ^a„,fcy'', a„,„ ^ 0. (5.13) 

fc=0 

Obviously ao,o = 1, a = ai,o and b = ai i. The orthogonality of the eigenfunctions {4>n} im- 
plies that {Pn[ri{x)^} are orthogonal polynomials in ri{x) with respect to the weight function 
00 (3;)^. Then we can show the following: 

%^^ = ^ + ^^r,H(x), r,^^\x)'='±vix + k), (5.14) 

0W(x) = 0H(x)pW(x) (n>.), (5.17) 

where 'Pn'(a;) is a symmetric polynomial of degree n — s in ri{x),rj{x + 1), . . . , ri{x + s) and 
satisfies the recurrence relation 

^ a.-M-.^-;'W-P^;''(x+ l)^ 

"^^ as^s ?7[^-il(x) - r7[^-i](x + 1) ^ ^ 
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Obviously Vn\x) = Pn{ri{x)) and vl^^x) = 1. From flS.lSp . Vn\x) is explicitly expressed as 

^ n s 

V^\x) = —Y,an,k J2 n^(^+^')'^ {n>s>0). (5.19) 

■^'^ k=s ko,...,ks>0 j=o 

For the shape invariant case, we have Vn\x) oc Pn-sivix] A + sS); A + sS). 
5.1.3 Casorati Determinants 

Here we prepare several formulas including various Casorati determinants and the sinusoidal 
coordinates. 

For any polynomial in rj, {Pn{ri)} {Pn{v) = c„77'^ + (lower degree terms)), a set of variables 
{f]j}, and a set of non-negative distinct integers {n^}, it is easy to show 

d^^iPnM)),<,,k<m= n iVk-V3)-RiVl,---,Vm), (5.20) 

l<j<k<m 

where R{rii, . . . , rjm) is a symmetric polynomial of degree YlT=i ~ \m{m — 1) in r^i, . . . , rjm- 
Especially, for (ni, . . . , n^) = (0, 1, . . . , £), the polynomial R becomes a constant and we have 

det(P,_i(r^,))^^^.,<,^^ = nc. ■ n (^fc-^.)- (5-21) 

fc=0 l<j<fc<^+l 

For all the sinusoidal coordinates rj{x) studied in [H |9], we have 

rj{x + a) + ?7(x — a) = (a polynomial of degree 1 in ?7(x)), (5.22) 
rj{x + a)rj{x — a) = (a polynomial of degree 2 in ^^(x)), (5.23) 

with G M. Hence a symmetric polynomial in T]{x),ri{x + l),...,r7(x + £) becomes a 
polynomial in 77 (x + |). Moreover the sinusoidal coordinate ri{x] X) satisfies 

r]{x + a; A) = (a polynomial of degree 1 in ri{x; A + 2ad)), (5-24) 

with G M. Therefore a symmetric polynomial in r]{x]X),r]{x + 1; X), . . . ,ri{x + A) 
becomes a polynomial in ?7(x; A + £6). The sinusoidal coordinate ?7(x; A) satisfies 

r](x + a; \) - r](x; X) , ^ , ^^^^ r.^\ 

\ \. = ^(x; A + (a - 1)(5). (5.25) 

r][a; A) 
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For the definition of the auxihary function ip{x] A), see eqs. (4.12) and (4.23) in [Tj. Let us 
define the £-th auxihary function (pi{x\ A) as 

ife{x;X)= [[ = II if{x + j;X+{k- J -1)5). [5.26) 

o<j<k<e-i 3t ) o<j<k<e~i 

Then, for any polynomials of degree n in ri{x; A), Pn{x; A) =^ Pn{ri{x; A); A), we have 



(5.27) 



W[Pn^,Pn^, . . . ,/'„^](x; A) ^ /a polynomial of degree YlT=i nk-\m{m-l 
V9m(x;A) \ in r7(x; A + (m - 1)5) 

and fl5.2ip implies 

£ Ik 

w[Po, A, . . . , Pi\{x- A) = n cfc(A) ■ n n ^) ■ (^.28) 

fc=o fc=ii=i 

Note that, for the orthogonal polynomials studied in [T], we have (lA.lGp . 
5.2 Formulas for modified Crum's theorem 

In this subsection we consider the simplification of the formulas in §|l]when the eigenfunctions 
consist of the orthogonal polynomials {Pn}'- 

<i)n{x) = M^)Pnivix)), Pnix) = Pn{vix)), (5.29) 

where Pn{v) is a polynomial of degree n in rj, satisfying the boundary condition 77(0) = 

dan. 

First we recall that the ground state (f)o{x) is annihilated by A, 



^fB{x) Mx) = ^/D{x + l)Mx + !)• (5-30) 
By using this and (13.281) . eqs. (I4.49p - (l4.50p become 

k=i k=i W[Pd,,...,PdA[x) 

Y\Dd, d,{x + i+l-k) = T]D{x + k)- W[P,„...,P,J(x) 

ti L\ W[P,„...,P.J(x + l) ^ ^ 

From these and f l4.45p -( l4.48p we obtain 



(x) = (-l)'Mx) 



f[B{xlk 1) W[P.„...,P.,,P„.](x) 
\ k=i JW[P,„ . . . , P,,](x)W[P,„ . . . , P,,](a; + 1) 



(5.33) 
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\ k=i VW[Pdi, . . . , PdJ(x)W[P,„ . . . , P,J(x + 1) 

B(x) = Bix + i) nP^^^---,P^M^) W[P,„...,P,„P,](x + l) 
W[P,,,...,PrfJ(x + l) W[Prf„...,P,„P^](x) ' 

Dix) -Dix) "^^^'^--- ^'^^^"^ + 1) W[P,„ ■ ■ ■ , P,,,P,](x - 1) 
W[P,,,...,P,J(a;) W[P,„...,P,„Pj(x) 

The eigenfunction 0„(x) of this P-deleted system has the following structure, 

T , , T , , W[Prf,,...,P,„P„](x) 

= (pJx) X — ^ - ' . • (5.37) 

^ ^ '^^ ^ W[P,„...,Prf„P^](x) ^ ^ 

In the following we assume that ri{x) satisfies fl5.22l) - fl5.25l) . Then the above formulas fl5.33p - 
fl5.36p are simplified thanks to fl5.27p : 

Pnfa(x)) , , def 



X- A) = ^x; A) X '^'')''Zii , Viix) = r]{x- X + id), (5.38) 



^ ' ^ ' V{r],.i{x)) V^{r],{x)) v,{x-\) <^,+i(x;A) ^ ' 

where the deformed polynomials V{ri£_i), VniVi) ^md V^{rji) are defined as (we suppress 
writing A-dependence explicitly) 

P(r7,_i(x)) = ¥.,(x; A)-i W[P,„ . . . ,P^J(x), (5.41) 

VM^)) = A)-i W[Prf„ . . . , P,„ P^](a;), (5.42) 

P„(r7,(x)) = <^,+i(x; A)-iW[P,,, . . . , P<,„ P„](a;), (5.43) 



and they are of degree \V\ - - I) , I'D] + - i£(£ + 1) and |P| + n - |£(£ + 1) in r^^.i, 
and ?7£, respectively. Here |P| is defined by \T>\ =^ X]j=i'^j- -^Y construction the deformed 
polynomials V{rii_i{x)) and V^{rii{x)), which appears in the denominator of fl5.38p and in 
(15. 39p - ( 15^^01) ■ are positive definite for x = 0, 1, . . . , Xmax- 

The ratio of the polynomials Vn{ve{x)) /V^{r]£{x)) are eigenfunctions of the similarity 
transformed Hamiltonian (a second order difference operator) "H: 

n = X)-' o n o 0^(s; A) = P(x)(l - e^) + D(x)(l - e"^), (5.44) 
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The orthogonality relations of {0n} (I4.5ip are now rewritten as the orthogonality relations 
of the polynomials {'P„(?7^(x))}, 



i'ix; XfVn{'ni{x))Vm{Vl{x)) = - £{dj)) ■ -;^^nm {n,m G {0, 1, . . . , ninax}\^^), 

(5.46) 



0^2 

x=0 j=l " 



in which the new weightfunction is defined by 

We have demonstrated the following: When the modified Crum's theorem is applied 
to the shape invariant (exactly solvable) Hamiltonian systems corresponding to the (g)- 
Askey scheme of hypergeometric orthogonal polynomials Pn{i]{x)) [Ij, the resulting system 
is no longer shape invariant but it is exactly solvable since all the eigenvalues and the 
corresponding eigenf unctions are explicitly given. The eigenfunctions consist of deformed 
polynomials Vn{f]i{x)) f l5.43p . They inherit the sinusoidal coordinate of the original system 
and they are orthogonal with respect to the weight function given as the square of the 
deformed groundstate eigenfunction ?/;(x; A)^ fl5.47p . These deformed orthogonal polynomials 
'Pn{ve{^))j although forming a complete set of eigenpolynomials, have vacancies in the degrees 
corresponding to the deleted levels V = {di,d2, ■ ■ ■ ,de}. Therefore they do not satisfy 
three term recurrence relations and they are not called orthogonal polynomials in the strict 
sense. When another (modified) Crum's theorem is applied to such a system, most of the 
simplification formulas (due to polynomial eigenfunctions) in this section are still valid with 
some adjustments. 

6 Dual ChristofFel Transformations 

In this section we will show that the simple case {i.e. /i = 0) of the modified Crum's the- 
orem applied to an orthogonal polynomial Pn{ri{x)) provides the well known Christoffel 
transformation [21], |3j for the corresponding dual orthogonal polynomial Qx{S{n)) fl2.2ip . 
Therefore we will call it the dual Christoffel transformation. Corresponding to the deletion 
of the levels or degrees n E V = {di, . . . ,de} of the polynomials Pn{ri{x)), the positions 
n E V = {di, . . . , de} are deleted from the dual polynomial Qx{£{n)). Here we assume that 
the set of deleted levels does not contain 0, the original ground state, i.e. ^ V. This also 
means that £ is even and the modified groundstate has also the label 0, /z = 0. The dual 
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Christoffel transformation needs the framework of the discrete quantum mechanics and it 
cannot be apphed to the general orthogonal polynomials. 

The deformed dual polynomials Qx{£) are defined by the three term recurrence relations 
in terms of B{x) flOSjl and D{x) f l536|) : 

B{x){Q,i£) - Q.+i(^)) + D{x){Q^i£) - = SQ.iS) (x = 0, 1, . . . , x^^x), 

Qo{£) = 1, Q_i(^) = 0. (6.1) 

One simple consequence of the above recurrence relation is the universal normalisation 

Q,(0) = 1 (x = 0,l,...,w). (6.2) 

Assume that ri{x) satisfies fl5.22p - fl5.25p . At the eigenvalues £ = £{n), the Schrodinger 
equation fl5.45p for the ratio of the deformed polynomials Vn{Vi{^)) fPoiVii^)) is identical 
with the above three term recurrence relation and the duality relation holds: 

in which 



PnQxiSin)), X e [0, 

•^max] ) ^ 



by the normalisation condition. We have verified the second equality for all the examples 
in [1], see flA.16p -f rA.17p . Thus Qx{£) are genuine orthogonal polynomials, satisfying the 



orthogonality relation dual to fl5.46p : 



nG{0,l,-,nmax}\X' i = l 



\2 



- ^2^, J ^xy (a;,y = 0, 1, . . . ,a;max)- (6.5) 

Note that 11^=1 (i^(?^) — £-{dj)) > due to (14. II) . It is obvious that the values of Qx{S{n)) at 
^positions' n eT) = {c?i, d2, ■ ■ ■ , di} do not enter the orthogonality relation or the Schrodinger 
equation fl5.45p . In other words, the values of Qx{^^{n)) at n E V are not defined but the 
degree x runs from to Xmax without any hole. This transformation of the dual orthogonal 
polynomials, Qx{^{i^)) ^ Qx{^{i^)) is the well known Christoffel transformation PT | [5| U7 \ 
. The transformation of the orthogonality weight function is read from (16. 5p : 

e 

dl^dll[{S{n)-£{d,)). (6.6) 

i=i 
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The inverse of the Christoffel transformation is called the Geronimus transformation [27], 
which adds new 'positions'. Its effects can be practically incorporated by the dual Christoffel 
transformations with redefinion of the parameters. 

The dual Christoffel transformation has the merits that the formulas determining B{x) 
and D{x) (I5.35P and (I5.36P are universal, concise and algorithmic compared with the orig- 
inal Christoffel transformation, which should be performed specifically for each case, the 
polynomials and the set of deletion V. While the dual Christoffel transformation requires 
the framework of the discrete quantum mechanics, the original Christoffel transformation 
is defined for general orthogonal polynomials, i. e., those having the three term recurrence 
relations only without the difference equation, nor the sinusoidal coordinate. The elementary 
Christoffel transformation is defined by 



for a parameter a G M, which is not constrained by the sinusoidal coordinate. Correspond- 
ing to the ^-deletion above, one considers multiple {i times) applications of the elementary 
Christoffel transformations. The normalisation weight function for the polynomials {P„(x)} 
is, in general, not known explicitly and it usually contains continuous measure. Thus de- 
termination of the positivity of the resulting weight function generically involves a moment 
problem. For the polynomials belonging to the discrete quantum mechanics PQ [S] , explicit 
forms of multiple Christoffel transformations are known for some specific cases. For example, 
see [28l EH] for the Racah polynomials and others. 

7 Summary and Comments 

Crum's theorem and its modification a la Adler are formulated for the discrete quantum 
mechanics with real shifts. They are slightly more complicated than those in the ordinary 
quantum mechanics or in the discrete quantum mechanics with pure imaginary shifts, partly 
because the size of the Hamiltonian matrix, or the range of the x variable (xmax), is reduced 
by one in each deletion. Another source of complications is the fact that in the generic cases 
the sinusoidal coordinate //(x) depends on the parameters A, which changes at each deletion. 
In the discrete quantum mechanics with real shifts, the modified Crum's theorem generates 
the Christoffel transformation of the dual orthogonal polynomials. Very special and simple 
examples, in which all the excited states from the first to the i-th are deleted (see Fig. 3), 



def Pn+l{x) - AnPn{x) 




(6.7) 



■n 



X — a 
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are presented explicitly in Appendix for more than two dozens of orthogonal polynomials 
discussed in [1]. 

One of the motivations of the present research is the connection with the infinitely many 
exceptional orthogonal polynomials [301 EH [321 |33l [M]. As explained in some detail in §4 of 
|19j . the insight obtained from the explicit examples of the application of modified Crum's 
theorem was instrumental for the discovery of the infinitely many exceptional Laguerre and 
Jacobi polynomials in the ordinary quantum mechanics [32] and the exceptional Wilson and 
Askey- Wilson polynomials in the discrete quantum mechanics with the pure imaginary shifts 
[M] . The explicit examples in the present paper are also very helpful for the discovery of the 
(infinitely) many exceptional orthogonal polynomials of a discrete variable, for example, the 
exceptional g-Racah polynomials, etc. They will be derived and discussed in detail elsewhere 
[35]. 
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A Special Examples 

In Appendix we present very special and simple examples of an application of Adler's theo- 
rem, in which the eigenstates 0i, 02, • • • , 0£ are deleted. Similar examples of the application 
of Adler's theorem in discrete quantum mechanics with pure imaginary shifts were given in 
Appendix B of [19] . In the ordinary quantum mechanics, some simple examples of the same 
sort were demonstrated in Appendix A of [19] and [36] . 

In this case, P = {di, d2, . . . , de} = {1,2,..., £}, that is, the modified groundstate level 
is the same as that of the original theory /i = 0. This means \V\ = + 1) and the 
polynomial "P^ fl5.42p is a constant. The degrees of the deformed polynomial V fl5.4ip and 
Vn (I5.43P are i and £ + n, and they are proportional to the deforming polynomial and the 
deformed polynomial P^^„ flA.24p . respectively. The duality (16. 3p is now a relation between 
two polynomials Vnive) and Qx{S). 

The main results are the unified expressions of the functions B(^{x) flA.29P and D(^{x) 



(lA.30p . which specify the Hamiltonian "H^ together with the eigenfunctions (p^^ (lA.3ip and 
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4>e,n{x) (]A.32p . The unified expressions of tlie deformed eigenpolynomials Pi^n{x) are given 
in (IA.43p . An important ingredient specific to eacfi type of polynomials is the deforming 
polynomial C,i{x), which is listed in IA.3I The situation is illustrated in Fig. 3, which should 
be compared with Fig. 2a depicting the generic case discussed in section HI 



e 

£-1 



2 
1 


original modified after i + 1 

system system steps of Crum 

Figure 3: Special case 




The black circles denote the energy levels, whereas the white circles denote deleted energy 
levels. We write "H = 1^12. ..i, (pn = <Pi2...£n, -A = Ai2...£, B = Bi2...e, D = Di2...e etc. as 
Tie, (pe^n, Ae, Be, etc. This Hamiltonian T-Li = AIAi is hermitian for even i but may 
be non-hermitian for odd £. Since algebraic formulas such as the Casoratians are valid for 
even and odd i, we present various formulas without restricting to the even £. The original 
systems are shape invariant but the . . . , 0^)-deleted systems Tie are not. The rightmost 
vertical line in Fig. 3 corresponds to the Hamiltonian system T-L'i = AeA\ = 7/^^^^', which is 
shape invariant and it is obtained from T-Li by one more step of Crum's method. 

We apply Adler's theorem to the shape invariant, therefore solvable, systems whose 
eigenfunctions are described by the orthogonal polynomials studied in [1]; i.e. (g)-Racah, 
(g)-(dual)-Hahn, etc. Hereafter we display the parameter dependence explicitly by A, which 
represents the set of the parameters, and we follow the notation of [Tj. This is to emphasise 
that the basic structure of Crum's theorem and its modification discussed in the main sections 
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are well-founded for the most generic systems. 



A.l The original systems 

Here we summarise various properties of the original Hamiltonian systems studied in p[| to 
be compared with the specially modified systems to be presented in IA.2[ All the examples 
are shape invariant and exactly solvable. Let us start with the Hamiltonians, Schrodinger 
equations and eigenfunctions (x = 0, 1, . . . , Xmax, a;max = or oo, n^a_^ = N or oo): 

^(A) = y/B{x; A) - e^y/D{x; A), ^(A)^ = y/B{x; A) - ^^(x; A) e'^ (A.l) 

n{X)=A{X)U{X), (A.2) 

?/(A)0„(x; A) = £:(n; A)0„(x; A) (n = 0, 1, . . . , n^ax), (A. 3) 

0„(x; A) = 0o(x; A)P„(r/(x; A); A) (P„(0; A) = 1, P„(x; A) = P„(r/(x; A); A)) . (A.4) 

The explicit forms of the set of parameters A, the potential functions B{x; X) and D{x] X), 
the energy eigenvalues £{n;X), the sinusoidal coordinate ?7(x;A), the ground state wave- 
functions (poix; A) and the eigenpolynomials Pn{f]{x\ A); A) are given in [I]. The groundstate 
wavefunction 0o(a;; A) is annihilated by ^(A), A{X)(I)q{x; A) = 0, and given by 



6o(x; A) 



The systems are shape invariant, 

A{X)A{X)^ = KAiX + 5)U(A + 6)+ £{1; A), (A.6) 

where the explicit forms of 6 and k are given in [I]. The action of .4(A) and .4(A) ^ on the 
eigenfunctions is 

.4(A)0„(x; A) = J-_ /„(A)0„_i (x; X + S), (A.7) 
a/-«(U; A) 

.4(A)V„_i(x;A + (5) = ^/B{0;X) 6„_i(A)0„(x; A), (A.8) 

where /n(A) and 6„_i(A) are the factors of the energy eigenvalue, S{n; A) = /„(A)6„_i(A), 
and their explicit forms are given in yj. The forward and backward shift operators J^{X) 
and B{X) are defined by 



J-(A) = y/BiO;X) 0o(x; A + 5)"^ o ^(A) o 00 (x; A) 
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B{0;XMx;X)-\l-e'), (A.9) 



= A) - D(a;; A)e-^)^(a;; A), (A.IO) 

and their action on the polynomials is 

T{X)Pn{v{x; A); A) = /„(A)P„_i(r/(x; A + 5); A + 5), (A.ll) 
S(A)P„_i(r7(x; X + d);X + 6) = 6„_i(A)P„(r/(a:; A); A). (A.12) 

For the definition of the auxiliary function ip{x] X) see eqs. (4.12) and (4.23) in [1]. Their 
explicit forms are also given in § 5 of [1] . 
The orthogonality relation is 

(f)o{x] XfPn{ri{x; A); X)P,n{'n{^\ A); A) = (n, m = 0, 1, . . . , n^ax), (A.13) 

where the explicit forms of the normalisation constants dn{X) are given in [1]. 

The coefficient of the leading term c„(A), which appears in Pn{y] X) = Cn{X)P^°^^'^{y] A), 
is given by 

11^=1^0; A) J-^A 5(0;A + j(5) 

for all the examples in [Ij. Especially n = 1 gives the relation 

P(0;A) _ -c,(A)-^ fAl^l 
^(1;A) - r,(l;A) ' ^^-''^ 

which expresses the important relation among the basic quantities 5(0; A), £^(1; A), //(I; A) 

and ci(A), see eq. (4.55) of [1]. The Casorati determinant ( 15.27^ becomes 

<^„(x;A)-iW[P„„...,P„„](x;A) 
rr' rx^ U,<,<,<JS{n,;X)-£{nr,X)) "^'^ ,^ ^ , 

= (-l)(^).-(T) n gK;A)-g(n,;A) _ 

V ; 11 B(0;X + (j -1)6) ("1'-'"-)^' ^' ^ ^ 

for all the examples in p. Here P(ni,...,n™)(a;; A) is a polynomial of degree XlfcLi ^fe — 1^(^— 1) 
in ri{x; A + (m — 1)<5) and satisfies the normalisation 

PK,...,n„)(0;A) = l. (A.17) 
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A. 2 The (01, . . . , 0£)-deleted systems 

For this very special case, the potential functions ( I5.35p -f l5.36p . the Hamiltonian f l4.39p - 
f l4.40p and the Schrodinger equation f l4.43p of the modified system are: 

B. i.; A) ^ Bi.; A) ^ Bi. + t, ^) w[A, . . . , P.](x + 1; A) W[A, • • • , P. Po](x; A) ' ^^'^^^ 

A) ^ A) ^ Dj.; A)^^^^ " " ' ^^^^ + V " " ' ^'^^^ " (A.19) 

^ ^ ^ ^ ^ ^ W[Pi,...,P,](x;A) W[Pi,...,P,,Po](x;A) ^ ^ 

^,(A) = ^(A) = ^/B,ix;X) - eVA(a:;A), 

^,(A)t = ^(A)t = v/5,(x;A) - ^DfXx;X) e~^, (A.20) 

•H,(A) = •H(A) = ^,(A)U,(A), xLx = N-i or oo, (A.21) 

?/,(A)0,,,(x; A) = f (n; A)0,,,(x; A) (n = 0, £ + 1, £ + 2, . . . , n^^). (A.22) 

Note that the Hamiltonian "Hi is an (x^^^^ + 1) x (^^max + 1) matrix He = {'He-x,y) ix,y = 
0, 1, . . . , x^j^x)- The explicit expression of De{x; A) can be analytically continued and gives 
^^(O; A) = 0. We have B^^x^^g^^) = for the finite case. From fl5.33p we have 

e-i 

(j)e,n{^; A) = ^n{x; A) = (-l)Vo(a:; A) JJ ^B{x + k;\) 

W[Pi, . . . , Pi]ix; A)W[Pi, . . . , P£](x + 1; A) 

Let us introduce the deforming polynomial ie{x; A) =^ ^e{rj{x; A + — 1)5); A), which is 
a polynomial of degree £ in //(x; A + (£ — 1)<5), and the deformed polynomial P£,„(x; A) =^ 
Pe,n{v{Xi A + M); A) (n = 0, £ + 1, £ + 2, . . .), which is a polynomial of degree n in r7(x; X + i6), 
in terms of (IA.16P as follows: 

ie{x; A) = P(i,2,...,^)(x; A), P£,„(x; A) =^ P(i,2,...,^,n)(x; A). (A.24) 

Then we have 

6(0; A) = 1, P,,„(0;A) = 1. (A.25) 

We set Pi^n{x; A) = for n = 1, . . . , For the g-Racah case, which is the most generic case, 
the deforming polynomial C,e{x; A) has the following form, 

g-Racah : ^e{x; A) = Pe{~x; t(A + (£ - 1)5)), t(A) = -A, (A.26) 
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and this is indeed a polynomial in ri{x; \ + {i — 1)6) because r]{—x; —A — {i — 1)6) = 
r]{x; X+{e-l)6){dq^-^)-\ For the other cases the deforming polynomials are obtained from 
this in certain limits. However these limits are not so trivial. So we present the explicit 
forms of C,i{x] A) in IA.3I Note that the deforming polynomial C,e is related to Pi^n as 

6(x;A) = P,-i/x;A), (A.27) 

and satisfies the recurrence relation: 

B{0; A + e6)^e+i{x; A) = B{x; A + e6)(p{x; A + i6)^i{x; A) 

- D{x; A + i6)ip{x - 1; A + i6)ii{x + 1; A). (A.28) 



By using these quantities, eqs. (lA.lSp . (1A.19P and (1A.23|) become 



B,{x- A) = ^'Bix; A + i6)jP^^. (A.29) 

OXx + l;X) 

De{x; A) = k'D{x; A + i6) ^'^^ + ^A.30) 

^iix; A) 

C{i,X) 0o(x;A + M) 

^5(0, A + M) ^^,(a;;A)e£(x + l;A) 

0,,„(x; A) = 0,,o(x; A)(-l)^ TT ^^""^ ^ ~ f.^^' ■ A), (A.32) 



where C (£, A) is defined by 



For even £, the deforming polynomial C,e{x; A) = ^e{ri{x; X+{i — l)6); X) is positive at integer 
points X = 0, 1, ... , x^a^ + 1. It has either no zero in the interval < x < x^^^^ + 1 or when it 
has zeros, an even number of zeros appear between two contiguous integers. The deformed 
polynomial Pt^nijJ] A) {n> ^ +1) has n — I zeros in the interval < y < Tji^xl^^^] X + i6) for 
even I. These we have verified by numerical calculation for lower ^ for all the examples in 
[1] . To the best of our knowledge, no general proof of this property has been reported. Note 
that the normalisation of does not affect H^. This system is not shape invariant. We have 
not chosen the normalisation like as 0o(O; A) = 1, namely 0£,o(O; A) = c(£^) ^ _^ 
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The operators ^^(A) and ^^(A)^ connect the modified system 'H^(A) = ^^(A)'''^^(A) to 
the shape invariant system n'^{X) = ^^(A)^^(A)^ = k.^+^U{\ + (£ + 1)8) + £{l + 1; A) = 
^/[^^^^(A), which is denoted by the rightmost vertical fine in Fig. 3. The n-th level {n > i+1) 
of the modified system Tie is iso-spectral with the — £ — 1-th level of the new shape invariant 
system 'H'f 

^,(A)0,,„(x; A) = /,,„(A)0„„,_i {x; A + (£ + , (A.34) 
A(A) Vn-f-i {x; A + (£ + 1)6) = be,n-i{\)(pi,n{^; A). (A.35) 

Here, fi,n{^) and 6£,„_i(A) are the factors of the energy eigenvalue, S{n; A) = /£,„(A)6£,„_i(A), 
and are defined by 

Note that £{n; A) = K^+^£{n -£-l;X + {£ + 1)6) + £{£ + 1; A) forn > £ + 1. The forward 
and backward shift operators J^^(A) and Bi{X) and the similarity transformed Hamiltonian 
T-Li^X), which act on the polynomial eigenf unctions, are defined by: 

J-,(A) = 00 (x; A + (£ + 1)6Y^ o A,{\) o 0,,o(a;; A) 

(l-e^), (A.37) 



v{x-\ + £6)0Xx + l-\) 

Bi{\) = 0^,o(x; A)-i o Ae{\)^ o 0o(x; A + (£ + 1)6) (A.38) 

= A + A) - D{x- A + + 1; A)e-^) <^(x; A + £5), 

H£(A) = 0^,o(x; A)-^ o H,(A) o 0,,o(a:; A) = i3f(A) J-,(A) (A.39) 

= k'{b{x- a + i6)-M^^{l - e^) + A + - e"^)) . 

Their action on the polynomials is 

J-,(A)A,„(x; A) = /,,„(A)P„_,_i(x; A + (£ + 1)(5), (A.40) 

-B,(A)P„_,_i(x; A + (£ + 1)5) = 6,,„_i(A)P,,„,(x; A), (A.41) 

?^,(A)P,,„(x;A) =^(n;A)P,,„(x;A) (n = 0, £ + 1, £ + 2, . . . , n^ax)- (A.42) 

For n > i + 1, the above formula ( 1A.38P provides a simple expression of the modified 
eigenpolynomial Pe^n in terms of and the original eigenpolynomial P^: 

«:-iC(£, A)6,,„_i(A)P,,„(a:;A) 
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= B{x; A + A)(/?(x; A + M)P„_£_i(x; A + + l)d) 

- D{x; A + id)ie{x + 1; A)^(x - 1; A + - 1; A + + 1)(5). (A.43) 

The orthogonality relation is (x^^x = — ^ or oo and rimax = or oo): 

0^,o(a;; A)^P^,„(x; A)P£,m(a;; A) = - — -— 6nm {n,m = 0,£+l,£ + 2, . . . ,nmax)- (A.44) 
^ a£,n(A)^ 

Here the normalisation constant d^^nW is 

which is a consequence of fl5.46p and positive for even £. The weight function 0£,o(a^; A)^ is 
positive definite for even i. 

A. 3 Explicit forms of the deforming polynomial ^£ 

As shown in the preceding subsection, the various quantities, the functions Bi{x) (lA.29p . 
Di{x) flA.SOp . which specify the Hamiltonian and the eigenf unctions (p^^ flA.311) . 0£,„(a;) 



(]A.32p are determined by the deforming polynomial ^^{x] X) = ^^(?7(x; A + (^— 1)5); A). This 
deforming polynomial ^^(x; A) is positive at integer points a; = 0, 1, ... , x^^^-^ for even i. Here 
we give the list of the explicit forms of the deforming polynomial C,e{x', A). The list contains 
the name of the polynomial {e.g. Racah) and the corresponding section number of [Tj together 
with that of Koekoek-Swarttouw's review [5], the definition of the twist operator t and the 
form of the deforming polynomial C,e{x] A) in terms of the original polynomial Pe{x] A) (the 
Racah polynomial): 

name (§ of [1], § of [5]): 

Racah (5.1.1, 1.2): t(A) =^ -A, 

Ux; A) = Pei-x; t(A + {£ - 1)6)), (A.46) 
Hahn (5.1.2, 1.5): t(A) = -A, 

6(x; A) = Pe{-x; t(A + {£ - 1)6)), (A.47) 
dual Hahn (5.1.3, 1.6): t(A) = -A, 

Ux; A) = Pe{-x; t(A + {£ - 1)6) + (0, 2, 0)), (A.48) 
Krawtchouk (5.1.4, 1.10): t(A) = {p, -N), 
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Ux; X) = Pt{-x-i{\+{i-l)6)), 
g-Racah (5.1.5, 3.2): t(A) = -A, 
Ux-\)^Pt{-x-i{\+{l-l)5)), 



def 



g-Hahn (5.1.6, 3.6): t(A) = -(6, a, N), 

Ux] \) = Pi{x-N + £- 1; t(A +{£- 1)6)) 



:-l)Vg^^(^-i)(6;g), 



dual g-Hahn (5.1.7, 3.7): 

'g~^, g^, a~^6~-^g~^+^" 



6(2;; A) 



3<^2 



a-ig-^+i,g^-^+i 



g;6g 



quantum g-Krawtchouk (5.1.8, 3.14): 

n „-x+N-i+i 



^^(a^; A) 

dual quantum g-Krawtchouk (5.1.8, — ): 



q;q){pq;q)i, 



^^{x;X)^2cpiy^^_,^^ 



—S r 

q , q 



q;pq''^), 

g-Krawtchouk (5.1.9, 3.15): t(A) = -A, 

6(x; X)^P,{x-N + e- 1; t(A + {£ - 1)6) + (-2, 0))(-l)VV, 
dual g-Krawtchouk (in the standard parametrization) (5.1.9, 3.17): 



6(2;; A) = 301 



q , c q 



-,N-e+i 



'■,cq^), 



affine g-Krawtchouk (5.1.10, 3.16): 



„-x+N-e+i 
y ) y 



^N-e+l 



i)ygl^(^+i) 



t{\)'^'-{b,a,N) 



^e{x; A) = 201 
alternative g-Hahn (5.3.1, 

^e{x;X)^P^{x-N + e-l;t{X+{e-l)6)) 

a {b;q)e 

alternative g-Krawtchouk (5.3.2, — ): t(A) = —A, 

6(x; X)^Pe{x-N + £-l;t{X+{£- 1)6) + (-2, 0)) (-1) V'?" 



(pq ; q)e 

-l)'q-¥('-'\a;q)e 



alternative affine g-Krawtchouk (5.3.3, — ): 

-x+N-e+i 



^e{x; X) = 201 



g %g 



q 



N-i+l 



q;pq- 



5+^+1 \ . 



^ ipq]q)i 



Meixner (5.2.1, 1.9): t(A) = (-^, c), 
6(a;;A) = P,(-a;;t(A+(£-l)5)), 
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Charlier (5.2.2, 1.12): t(A) = -A, 

little g-Jacobi (5.2.3, 3.12): 1(A) = -A, 

Ux; A) = P,{x + b' + i; t(A + (i - 1)6) - (2, 2))a-%-'q-'^'+'\ b = q'' , 

dual little g-Jacobi (5.2.3, — ): 

- /q-^,q'',a-^b-\-''-^ 

g-Meixner (5.2.4, 3.13): 



ie{x; A) 



2Vl 



—f T 

q , g 

6-ig-^ 



g; -6 ^g^ ), 



little g-Laguerre/Wall (5.2.5, 3.20): 



ii{x; A) = 1^ 



a ^g" 



q;a-Y)i-l)'a-Y^A^^'\aq;q)e, 



Al-Salam-Carlitz II (5.2.6, 3.25): 



iiix; A) 



2(pi 



— f T 

q , g 




q;a ^q^ ^' 1 , 



alternative g-Charlier (5.2.7, 3.22): 



6 (a;; A) 



g , —a g~ 



g; -ag 



x+2t 



(-a)-V, 



dual alternative g-Charlier (5.2.7, 



A) 



— / T —1 1— r— / 

g , g , —a g 



0,0 



g;g 



g-Charlier (5.2.8, 3.23): 



6 (a;; A) 



2<P0 



— / T 

g , g 



(A.62) 
(A.63) 

(A.64) 

(A.65) 

(A.66) 

(A.67) 

(A.68) 

(A.69) 

(A.70) 



The twist operator t is not listed when it is not used for the definition of ^g. 

A. 4 Supplementary data on dual orthogonal polynomials 

Here we present the extra data for various dual orthogonal polynomials which were not listed 
in [1]. They are the dual quantum g-Krawtchouk in § 5.1.8, the dual little g-Jacobi in § 5.2.3, 
the dual alternative g-Charlier in §5.2.7 in [1]. 
• dual quantum g-Krawtchouk: 



g^ = (p,g^), 5 = (0,-1), K = q-\ p > g-^ 
£(n;A) = g-"-l, r/(x; A) = 1 - g^ 



(A.71) 
(A.72) 
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Pn{rj{x;X);\) = 2(j)i^ 



q n^q X 



q-pq 



dn{Xf = 



N 



iQ;q)x{q;q)N-x (p ""-.q). 

{q\q)N p-nqn(n-l-N) 



{q;q)n{q;q)N-n {p ^q ^]q)n 
Ri{z;X)^{q--2 -q'^fz', z'^^'z + 1, 



(^(x;A) = g^ /„(A) = r"-l, 6„(A) = 1. 

dual little g-Jacobi: 

q^ ^ {a,b), 5 = (0,1), K^q, < a < g"\ b<q~^, 
£{n- A) = 1 - g^ 77(x; A) = (g-^ - 1)(1 - a6g-+^). 



Pn(?7(a;; A); A) = 30i 
^Q{x;Xf = 



g "jg "^^abq""^^ 
bq 

.2 



g;a ^g" 



2 (6g, abq ; g)^; a^g^"" 1 — abq^^^^ 



{q;q)n 



{q,aq;q)^ 

{aqT X 



1 — abq ' 



^ (;^)2 _ (&9;9)n^„„^„ {aq;q)c 



R,{z-X) = {q-i-q^fz\ z' = z-l, 
i^o(^;A) = (g-^-g^)V^ 

i?_i(z;A) = (g-5 -g^)2((l + a6g)z'2 + (l + a)z'), 
g~^ — abq^'^'^ 



{abq^; g)oo' 

def 



(^(x; A) 



1 — abq^ 



-, /„(A) = l-g", 6„(A) = 1. 



dual alternative g-Charlier: 

q'^ = a, S = 1, n = q, a > 0, 

£(n; A) = 1 - g", rj(x; A) = (g-^ - 1)(1 + aq^, 



PM^;Xy,X) = q''\<P, 



q ",q ^ 




q--a q ^ = g^of 



g;-a g 



0o(2^; A)^ 



^x^lx{Zx ]^ _|_ 



ag 



,2x 



dn{X) 



2 a-g2 



(g;g)^ 1 + a ' 

n^in(n+l) 



X 



(g;g)„ (-ag;g)oo' 



A.73) 
A.74) 

A.75) 

A.76) 
A.77) 
A.78) 
A.79) 

A.80) 
A.81) 

A.82) 
A.83) 

A.84) 

A.85) 
A.86) 
A.87) 

A.88) 

A.89) 
A.90) 



A.91) 
A.92) 

A.93) 
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R^{z;X) = (q--^ -q^^fz', z' = z - 1, (A.94) 

Roiz;X) = {q--^-q^)h'^, (A.95) 

R.,{z; A) = (g-^ - g5)2((i _ a)z'' + z'), (A.96) 

y(x;A) = ^ - + ag-+^ ^ /„(A) = i-g", 6„(A) = 1. (A.97) 

We note that the httle g-Jacobi polynomial eq. (5.193) and the alternative g-Charlier poly- 
nomial eq. (5.248) in [1] can be rewritten as 

little g-Jacobi : Pn(r/(a;; A); A) = 3^1 (^^ ^ ^ ^ . ^-i^^ j ^ 

f q~^ g~*^ — (ig'^ \ 
alternative g-Charlier : P„(r7(x; A); A) = 3^0 ( ' ' ; — "'^Q'^ ) • (A.99) 
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